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In this paper, we study a two-lane totally asymmetric simple exclusion process (TASEP) 
coupled with random attachment and detachment of particles (Langmuir kinetics) in 
both lanes under open boundary conditions. Our model can describe the directed motion 
of molecular motors, attachment and detachment of motors, and free inter-lane transition 
of motors between filaments. In this paper, we focus on some finite-size effects of the 
system because normally the sizes of most real systems are finite and small (e.g., size < 
10, 000). A special finite-size effect of the two-lane system has been observed, which is that 
the density wall moves left first and then move towards the right with the increase of the 
lane-changing rate. We called it the jumping effect. We find that increasing attachment 
and detachment rates will weaken the jumping effect. We also confirmed that when the 
size of the two-lane system is large enough, the jumping effect disappears, and the two- 
lane system has a similar density profile to a single-lanc TASEP coupled with Langmuir 
kinetics. Increasing lane-changing rates has little effect on density and current after the 
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density reaches maximum. Also, lane-changing rate has no effect on density profiles of a 
two-lane TASEP coupled with Langmuir kinetics at a large attachment/detachment rate 
and/or a large system size. Mean-field approximation is presented and it agrees with our 
Monte Carlo simulations. 

Keywords: Langmuir kinetics; TASEP; two lane; density profiles. 

1. Introduction 

Recently, driven diffusive systems have attracted the intere sts of physi cists because 
the systems show surprisingly rich and complex behavior !^^ 4 |5J6| ^ g a g^p^ 
model of driven diffusive systems, asymmetric sim ple exclusion processes (ASEPs) 
have been widely studied in chemistry and physics!-^ Moreover, ASEPs hav e been 
applied successfully in biology such as gel electrophoresis^ protein synthesis I 
mRNA tr anslat ion phenomena}!^ motion of molecular motors along the cytoskeletal 



A totally asymmetric simple exclusion process (TASEP) is regarded as the min- 
imal model of ASEPs, in which particles move along one direction. TASEP has 
received much attention in model ing intra-cellular molecular transport and related 
traffic P roblem3^U- 7 l 18 l 19 l 2U l 21 l ^ l in recent years. 

MacDonald et alff^l initially used a TASEP to simulate the dynamics of ribo- 
somes moving along messenger RNA chains. Kruse and SekimotoISl proposed a 
two-headed TASEP model to describe molecular motor traffic. In their model, each 
motor consists of two heads and each head is attached to a filament. Each filament 
is assumed to be moveable. The investigation indicates that the average relative 
velocity of filaments is a non-monotonic function of the concentration of motors. 
The density and current profiles of motors with different types of microtubule tracks 
(e.g., cylindrical geometry and radial geometry) have been investigated by Lipowsky 



An extension of single-lane TASEP, incorporating Langmuir kinetics (LK, the 
at tac hment and detachment of particles), has been presented by Parmeggiani et 
gjdni Their model is also referred to as the PFF model. Numerical results of the 
PFF model show that there are unexpected stationary regimes for large but finite 
systems. Such regimes are represented by phase coexistence in both low and high 
density regions, which are separated by a domain wall. 

A mod el, incorporating a TASEP, Lang muir kinetics and Brownian ratchet 
mechanismj 2 -^ is proposed by Nishinari et all^U to mimic the movement of the 
single-headed kinesin motor, KIF1A. A novel feature in their model is that there 
are three states (strong binding, weak binding and no binding) of a KIF1A, com- 
pared to two states (binding or unbinding) in previous models. Their model can 
capture explicitly the effects of adenosine triphosphate (ATP) hydrolysis as well 
as the ratchet mechanism which drives individual motors. The experimentally ob- 
served single molecular properties in the low-density limit are reproduced and a 
phase diagram is presented. 



filamentsH^LSI and the depolymerization of microtubules by special enzyme: 
well as vehicular traffic 




as 
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Most previous work on modeling molecular motor traffic deals with single-lane 
systems where particles can move forward or backward, or attach and detach to a 
bulk (a collection of sites in a lane except boundaries) . Obviously, the description 
of the traffic of motors would be more realistic if multi-lane asymmetric exclusion 
processes can be considered. Experimental observations^! have found that motor 
protein kinesins can move along parallel protofilaments of microtubules and they 
can jump between these protofilaments without restraint. 

Recently, Pronina and KolomeiskjESl proposed a two-lane traffic model to simu- 
late a TASEP with symmetric lane-changing rules between two lanes, but without 
Langmuir kinetics. The computational results suggest that lane-changing rates have 
a strong effect on the steady-state properties of the system. In particular, the par- 
ticle current of each lane will decrease and particle densities will increase with the 
increase of particle coupling. Pronina and Kolomeisky the n ex tended their model to 
a more general case where asymmetric coupling is applied!^ It is found that asym- 
metric coupling between lanes leads to a very complex phase diagram, quite different 
from symmetric coupling. There are seven phases in the TASEP with asymmetric 
lane-changing rates, in contrast to three phases found in the system with symmet- 
ric coupling. In addition, a new maximal-current phase with a domain wall in the 
intermediate coupling is reported in RefP 

Effects of synchronization of kinks (i.e., domain walls) in a two- lane TASEP 
without Langmuir kinetics is investigated by Mitsudo and Hayakawa!^ The asym- 
metric lane-changing rate between lanes is used. Moreover, different injection and 
ejection rates of particles at the boundaries of two lanes are also considered. The 
positions of kinks are reported to be synchronized, though the number of particles 
may be different on these two lanes. 

More recently, Jiang et alP^ introduced Langmuir kinetics into one lane of 
a two-lane system. This has shown that synchronization of shocks on both lanes 
occurs when the lane-changing rate exceeds a threshold. A boundary layer is also 
observed at the left boundary as a finite-size effect. In their model, attachment and 
detachment of particles are assumed to occur only on one of two lanes, not both, 
which is unable to realistically describe real two-lane or multi-lane molecular motor 
traffic. 

This paper will investigate the collective effect of attaching and detaching par- 
ticles on both lanes of a two-lane system with symmetric inter-lane coupling. In 
particular, we will focus on the finite-size effects of the system. The model de- 
scribed in this paper is directly motivated by the dynamics of molecular motors, 
for instance, unidirectional motion of molecular motors along fil ame nts}^ random 
motor (e.g., kinesin) attachments and detachments to filaments}^ and molecular 
motors freely changing to adjacent filaments!^ It is expected that the incorporated 
process of a TASEP, Langmuir kinetics and lane-changing can shed light on the 
study of the traffic of molecular motors and other particle traffic in biology as well 
as in other areas such as vehicular traffic. 

We investigate the effects of the following parameters on density and current 
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profiles: (i) Different lane-changing rates u>. We denote O = uN, where f2 is used 
to represent the number of lane-changing particles and N is the length of each lane 
(i.e., the number of cells in each lane); (ii) Different combinations of attachment 
(Qa = waN) and detachment (Qd = ^dN) rates (for simplicity, we assume that 
attachment rate £Ia is fixed but detachment rate flo varies); and (iii) Different sys- 
tem size N, i.e., the length of each lane. Note that as we only consider symmetric 
lane changing, these two lanes are equivalent. We compare results obtained from our 
model when ft > and = (where the two-lane system becomes two separated 
TASEPs coupled with Langmuir kinetics). In this way, we can more clearly distin- 
guish the effect of different values of £1 on a two-lane system. Finally, mean-field 
approximation (MFA) is presented and used to verify our Monte Carlo simulations 
(MCS). 

The paper is organized as follows. In Section II, we give a description of our 
two-lane TASEP model, considering attachment and detachment of particles in 
both lanes. In Section III, we present and discuss the results of our Monte Carlo 
simulations. In Section IV, mean-field approximation is presented and the quanti- 
tative agreement between the Monte Carlo simulations and the mean-field theory 
is confirmed. Finally, we give our conclusions in Section V. 

2. Model 

Our model is defined in a two-lane lattice of N x 2 sites, where N is the length of 
a lane. We assume that all particles move from the left to the right, as shown in 
Fig. 1 . Sites i = 1 and i = N define the left and right boundaries respectively, while 
a set of sites i = 2, N — 1 is referred to as a bulk. Particles in the system arc 
involved in the following processes: lane selection, particle injection into the first 
site of a lane, particle detachment from a bulk, particle movement along a bulk, 
particle lane-changing between bulks, particle attachment to a bulk, and particle 
ejection from the last site of a lane. Thus, we assume that a particle that is injected 
to a lane and moves along the lane may be involved in the following processes: 

• Lane selection: A lane is randomly chosen. 

• Particle injection: If the first site of the selected lane is unoccupied, a par- 
ticle will be injected to the site with probability a. 

• Particle detachment: A particle can leave the system with probability u)d- 

• Particle movement: If a particle cannot leave the system, it can move for- 
ward to the next site in the same lane provided that the next site is empty. 
If the next site is not empty, the particle cannot advance. In this case, the 
particle will check if it can change to the other lane. 

• Particle lane-changing: The particle may change to the corresponding site 
on the other lane with probability ui if the corresponding site on the other 
lane is unoccupied. 

• Particle attachment: If a site in a bulk is empty, a particle can attach to 
the site with probability ua- 
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• Particle ejection: If a particle reaches the last site of a lane, it will be ejected 
from the system with probability (3. 



a 
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lane 1 



a 
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o 



p 



lane 2 



to. 



to. 



Fig. 1. The schematic representation of a two-lane TASEP with symmetric inter-lane coupling, 
and the attachment and detachment of particles on both lanes. 



The above processes can be described using the following updating rules. As the two 
lanes are equivalent, they have the same updating rules. We only list the updating 
rules for lane 2. An occupation variable (r^j) denotes the state of the ith site in the 
£th lane, where tu = 1 (or t?^ = 0) corresponds to whether the site is occupied (or 
empty). Thus, ti^ refers to the state of site i on lane 2. The rules are: 

• Case i = 1. (i) If T2 ; i = 0, a particle enters the system with probability a; 
(ii) If T2,i = 1 and T2 : 2 — 0, then the particle in the site of (2, 1) moves into 
site (2, 2); (hi) If r 21 = 1 and r 2 ,2 = 1, then the particle in the site of (2, 1) 
stays there. No lane change occursQ 

• Case i = iV. If T2,n = L the particle leaves the system with probability @. 
No lane change occurs[j 

• Case 1 < i < N. (i) If t^^ = 1, the particle may leave the system with 
probability lu^; If it cannot leave the system, it moves into site (2,i + 1) 

a The simulation results show that they are essentially the same if lane-changing is allowed on the 
first site. 

b The simulation results also show that there are no differences on density profiles whether lane- 
changing is considered in the last site. 
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provided T2S+1 = 0. If the particle cannot advance, it may change to lane 
1 with probability u> if — 0; (ii) If Ti^i = 0, a particle enters the system 
with probability lo a- 

These updating rules for both lanes arc illustrated in Fig. 1. For simplicity, we use 
a symmetric lane-changing rule. It will be the next step of our work to investigate the 
effects of asymmetric lane-changing rates on the dynamics of interacting particles. 

Normally, there are three types of updating schemes!^ (i) l n parallel. Updating 
rules are synchronously applied to all sites. A typical application is in vehicular 
traffic flow; (ii) In ordered sequence. Positions of particles are updated in an or- 
dered sequential manner, e.g., from one end to the other end of a one-dimensional 
system; (iii) In random. A site is randomly chosen and it is updated according 
to updating rules. Here we use a randomly sequential updating scheme, i.e., the 

third scheme, which has b een widely used in the simulations of molecular motor 
traffic llHll»l20l21l25l2f)l28l 
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Fig. 2. (Color online) (a) and (b) show average density p(x)\ (c) and (d) show current J(x) for 
different values of Q. The system parameters are set to: TV = 1000, Qa = 0.3, (2^ = 0.1, a = 0.2 
and = 0.6. (b) and (d) are the locally enlarged figures of (a) and (c) respectively. 



3. Monte Carlo Simulations 

In this section, the results of Monte Carlo 
current (Jj) is given by: Ji = pi(l - pi+\). 
a = 0.2, (3 = 0.6, VL D = 0.2, K = 3 and 



simulations are presented. The average 
The parameters are set to N = 1,000, 
Qa=K*VL y). In simulations, stationary 
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profiles are obtained by averaging 10 5 samples at each site. The sampling time 
interval is ION. The first 10 5 iV time steps are discarded to let the transient time 
out. For the special case when fi = 0, the two-lane process is reduced to two 
separated processes: two one-lane TASEP with Langmuir kinetics. Features of a 
one-lane TASEP with Langmuir kinetics have been well studied in RefP 

Fig. 2 shows average density p(x) and current J(x) of one lane of a two-lane 
system for different values of fi. For the special case O = 0, the two-lane system 
reduces to two separated single-lane systems. When fi increases from to 1, we 
may expect the domain wall always moves in one direction (either to the left or 
right). However, an unexpected phenomenon has been observed in our simulations: 
the domain wall first moves to the left slightly, and then moves towards the right 
with the increase of fi (see Figs. 2(a) and (b)). This process is more clearly captured 
in Figs. 3(a) and (b), where the domain walls for fi < 0.4 are shown. One can see 
that the domain wall moves to the left first when fi increases from to 0.2, and 
then moves towards the right when fi > 0.2. We refer to this feature as the jumping 
effect to describe the domain wall (as well as current curve) jump to the left and 
then to the right caused by the particles jumping between the two lanes. We also 
find that there is a critical value, fi c w 0.2, to distinguish the movement direction 
of the domain wall under certain conditions. A similar effect occurs in RefP 




j i i i i i i i i i i_ 
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Fig. 3. (Color online) Average density p(x) for small values of Q. The system parameters are: 
N = 1000, fl A = 0.3, n D = O.l, a = 0.2 and /3 = 0.6. 



When Qa and fip increase proportionally, the jumping effect becomes weaker 
(comparing Figs. 2(a) and (b)with Figs. 4(a) and (b)). This implies that the jumping 
effect is also influenced by attachment and detachment rates. 

The jumping effect is also observed in a larger system, e.g., N = 10,000, (see 
Fig. 5) when fi is small. It can be seen that the domain wall moves left when 
fi < 0.3. When fi > 0.3, the domain wall moves right. For instance, there is an 
obvious movement towards right when fi = 0.4 (see Fig. 5(b)). The critical value 
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Fig. 4. (Color online) (a) Average density p(x); Current J(x) for different values of Q. The system 
parameters are: TV = 1,000, Qa = 0-6, £Id = 0.2, a = 0.2 and /3 = 0.6. (b) and (d) are locally 
enlarged figures of (a) and (c) respectively. 



fi c of Q varies with a system size, i.e., fi c changes from 0.2 to 0.3 when the system 
size increases from 1,000 to 10,000. 

We have also confirmed that the jumping effect in a two-lane homogeneous 
TASEP coupled with Langmuir kinetics is a kind of finite-size effect. When the 
system size is very large, e.g., N — 100, 000, in our simulations, the jumping effect 
almost disappears. 

We argue that we need to study and examine finite-size effects, like the jumping 
effect because the size of a real system is normally not very large. For example, 
kinesin protein motors are responsible for long-distance transport in a cell. The 
length that a kinesin protein motor travels from its origin to its destination is 
normally about 100 successive steps on microtubules and the step size of kinesins is 
about 8nmf2SlThis indicates that a system size below 1,000 is enough to realistically 
describe the movements of molecular motors. Even for a large system, a size N < 
10, 000 is normally adequate for simulation. Thus, a conclusion obtained based on 
the assumption of an infinite size may not be applicable to any real system. 

Comparing Figs. 2(a) and (b) with Figs. 4(a) and (b), it can be seen that the 
slopes of domain walls increase with the increase of Qa arid Qd proportionally. 
This is because when the difference between Qa and fin enlarges, more particles 
will attach to these two lanes rather than detach from them. 

Increasing the value of fl means more particles have chances to change to the 
other lane if they cannot move forward along the current lane. Upon further in- 
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creasing f2, it is found that the density decreases, and the locations of shocks shift 
right slightly (see Figs. 2(a) and 3(a)). 

When x > x^ax ( x max denotes the position where the density reaches maxi- 
mum), the differences between density profiles are not obvious (see Figs. 2 and 4). 
This suggests that lane-changing rate does not have any obvious effect on the 
system properties (e.g., average density and current profiles) after the maximum 
density is reached. 

In Figs. 2(c) and (d), and 4(c) and (d), current profiles are illustrated. It is 
found that (i) The current decreases with the increase of 17 until it approaches 
the maximum current (Jmax)- When pi+i is approximately equal to pi, the current 
Ji w pi(l - pi). Thus, the maximum current J m ax ~ 0.25. We denote by x™£ x the 
position where the current reaches the maximum; (ii) When x^ x < x < xf^" x i tne 
current increases when lane-changing rate f2 increases; (iii) When x > xf%™ x , the 
current profiles are almost the same for different f2. 




Fig. 5. (Color online) Average density p(x) for small values of Q in a large system. The system 
parameters are: N = 10000, Q A = 0.3, n D = 0.1, a = 0.2 and (3 = 0.6. 



We next investigate the effect of the values of K (K = Ha/^d) on average 
density and current profiles. We consider the situation in which £Ia is fixed but 
fl o varies. Note that the increase of K means the decrease of VI d- It can be seen 
that average density increases with the increase of the value of K (see Fig. 6(a)). 
Moreover, the shock moves left and the amplitude increases. In other words, the 
width of the transition region decreases. This may be explained as follows: when 
the value of K increases, the value of fljj decreases proportionally, so that the 
opportunities for particles to detach from the bulk become smaller. In this case, 
more particles will remain in the bulk. This leads to an increase of average density. 

In Fig. 6(b), the current first increases with the increase of the value of K till it 
reaches the maximum value, and then the current decreases with the increase of the 
value of K. In other words, the position with the maximum current moves left and 
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Fig. 6. (Color online) Average density p(x) for different values of K. The system parameters are: 
N = 1000, fl A = 0.3, n D = O.l, a = 0.2 and = 0.6. 



the amplitude increases (see Fig. 6(b)) upon increasing K. The current curves can 
be obtained from equation Jj = pi(X - Pi+\), it thus can be seen that the maximum 
current J ma x =0.25 when densities of two adjacent sites are equal to 0.5. When 
density is higher or lower than 0.5, the corresponding current is lower than 0.25. 




Fig. 7. (Color online) (a) Average density p(x); (b) current J(x) for different values of the system 
size. The system parameters are: Q = 0.2, Qa = 0.6, flo = 0.2, a = 0.2 and f3 = 0.6. 



Fig. 7 shows the density and current profiles for different system sizes. It is found 
that the width of the transition region decreases as the number of sites is increased 
(see Fig. 7(a)). Our simulation suggests the boundary layer shows a finite-size effect 
and will disappear in the limit of an i nfinite system. The finite-size effect shown in 
Fig. 7(a) is consistent with that in Refs.^EHln Fig. 7(b) one observes the amplitude 
in the left boundary enhances when the system size increases. 
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4. Mean-field Approximation 



In this section, a mean- field theory is developed. The occupation variable (rej) 
defined in section II is also used. The corresponding equation for the evolution of 
particle densities (r^j) in a bulk (i.e., 1 < i < N) can be written as 

fell = (T M _!(1 - T M )> - (T M (1 - T M+ l)> 
+ w(T2,iT2,i+l(l - T M )) - w(ri 1 iTi ) i + l(l - T 2 ,i)) 

+wa(1 -n,i) - ^(ti^), (l) 



^ = (r^i-xCl - 75,0) - (7*,i(l - r 2>i+1 )> 
+ w(ri,<ri,t+i(l - T2,i)) - w(r 2l iT2,i+i(l - n,;)) 

+Wa(1-T2,») -W/j(T 2 ,j), (2) 

where (•••} denotes a statistical average. The term w(ri i iTi j , + i(l — T2,»)) is the 
average current from site i in lane 1 to site i in lane 2, while k>(T2,jT2,i+i(l — 
is the average current from site i in lane 2 to site i in lane 1. The first two terms in 
Eqs. (1) and (2) correspond to particle movement; the middle two terms correspond 
to particle lane changing, and the last two terms correspond to attachment and 
detachment of particles. At the boundaries, the densities evolve as 

d(n t x) 



dt 

d(n tN ) 



= a(l-Ti, 1 )-(T t , 1 (l-Tl,2)), (3) 

= (n,N-i{l -tim)) -/3(n tN ), (4) 



dt 

where I can be 1 or 2. It can be seen that Eqs. (1) and (2) are equivalent when we 
use the same lane-changing rate between two lanes. Therefore, we only focus on the 
mean- field approximation in lane 2 (i.e., Eq. (2)). When N is large, i.e., N 1, 
we can transfer the coarse-grained stationary-state Eq. (2) to continuum mean- 
field approximation. First, we factorize correlation functions by replacing (r^j) and 
(r^i+i) with pi t i and pe t i+i, then we set 

Pl , l±1 = p(x) ± -- + —— + (_). ( 5 ) 
Substituting (J5]) into Eq. (2), we obtain 

^ = (2 P2 - 1)^ + q a (i - pi) - n DP2 
-n(i- Pl )p 2 2 + npl(i-p 2 ), (6) 

where t' = t/N, — coN, fl^ = io^N and fljj = ui^N. The term (ti <ti T2 <)) 

is replaced by pi,ipi,i+i{l — p2,i), and pi.i+i is approximated as pi^; similarly, the 
term (t2 iT2,i+i(l— t^)) is replaced by P2,ip2,i+i(l—pi,i), andp 2 ,i+i is approximated 



9, 2008 3:45 WSPC/INSTRUCTION FILE wangrl 



12 R. Wang et al, 

as p2,i- The boundary conditions become pt(x = 0) = a and pe(x = 1) = (1 — /?). As 
lane-changing is symmetric and two lanes are homogeneous, we have pi = p± = p 
(our Monte Carlo simulations show that densities in lanes 1 and 2 are the same). 
Thus, Eq. (6) reduces to 

^ = (2p-i)^ + n A (i- P )-n DP . (?) 

In the limit of t — » oo, the system reaches a stationary state with ^ = 0, Eq. 
(7) simplifies into 

(2p-l)^+(l A (l-p)-Sl D p = 0. (8) 

Eq. ([5]) is a first order differential equation and has been solved in RefP^ The 
low-density profiles in the bulk can be obtained by integrating the equation from 
the left boundary (p(0) = a) to a density pf. 

1 , n , s K-l, ,K-(K + l)pi n 

-[2(p,-a) + — — In |— 4— -fi|]. (9) 



f2 D (K+l) K + l 'K-(K + l)a 

Similarly, we can also integrate the equation from the right boundary (p(l) = 1-/3) 
to a density p r : 



Q D {K + l) iK r K + l 1 K-{K + l)p r 

where K = Ha/^d- The overall density profile across the system can be described 
by determining a shock position in these two profiles. In Fig. 8, we compare the 
results of mean-field approximation (MFA) with the results of Monte Carlo simula- 
tions (MCS). We find that the results of Monte Carlo simulations are in agreement 
with that of mean-field approximation. We note that nume rica l results of MCS in 
Figs. 8(a) and (b) are almost the same as Figs. 2 and 5 in Ref.pSlwhich also suggests 
that lane-changing rate £1 has no effect on density profiles of a two-lane TASEP cou- 
pled with Langmuir kinetics at a large attachment/detachment rate and/or a large 
system size. 



5. Conclusion 

In this paper, two-lane totally asymmetric exclusion processes coupled with Lang- 
muir kinetics on both lanes are studied using Monte Carlo simulations (MCS) and 
mean- field approximation (MFA). The results of Monte Carlo simulations agree 
with the results of mean-field approximation. The model is directly inspired by 
the experimentally observed movement of molecular motors which include motor 
advancing along filaments, random attachment and detachment, and free jumping 
between filaments. 

The system has mainly demonstrated the following complex behavior on two 
lanes. 
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1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 
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Fig. 8. Density profiles from Monte Carlo simulations (black line) and mean-field approximation 
(red dashed line), (a) a = 0.2, (3 = 0.1, k = 1,Q D = 0.1 and fl = 0.2; (b) a = 0.2, /3 = 0.6, k = 
3,tt D = 0.1 and U = 2; (c) a = 0.2, (3 = 0.6, k = 3,tt D = 0.2 and U = 100. The system size 
N = 1,000. 



• The jumping effect is observed, that is, the domain wall first moves left 
slightly, and then move towards the right with the increase of the lane- 
changing rate. This effect is a finite-size effect as it is not observed when 
the system is large, e.g., N = 100,000. On the other hand, increasing at- 
tachment and detachment rates proportionally will lead to the jumping 
effect becoming weaker. 

• After densities reach maximum, the increase of the lane-changing rate has 
little effect on the system properties such as average density and current. 

• When (attachment rate) is fixed, it is found that average density and 
current increase and then the current lowers upon decreasing the value of 
fin (detachment rate). 

• When system size iV is increased, the amplitude of the shock increases and 
average density increases. 

• Lane-changing rate has almost no effect on density profiles in a two-lane 
TASEP coupled with Langmuir kinetics at a large attachment /detachment 
rate and/or a large system size. 
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